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Abstract 

The top-down approach of engineering software integration is consid- 
ered in this parer. A set of advantages of this approach are presented, by 
examples. All examples are supplied by open source code. 



1 Introduction 

Historically engineering software integration is rather chaotic then planned. Dif- 
ferent CAD products had been integrated with CAE and CAM products, Mat- 
lab had been integrated with Simulink etc. This integration history could be 
schematically presented in figure [TJ 

However chaotic way is not optimal way. For example CH — I- programming lan- 
guage had been designed chaotically. Languages of next generation (Java, C#) 
had been designed by planned way. C++ does not support automatic garbage 
collection. But any present day project should have it. So big C++ projects 
contain smart pointers those provide automatic garbage collection. Since smart 
pointer is not intrinsic C++ feature there exists a lot of versions of smart point- 
ers. Integration of set of projets with different versions of smart pointers has a 
set of disadvantages. In this case integrated product is not clear, volume of its 
code is overloaded etc. Present day programming languages (Java, C#) sup- 
port Reflection but C++ does not support it. However Reflection is very useful 
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Figure 1: Chaotic integration of software products 



feature for any big project. So big C++ projects have different versions of Re- 
flection emulation. Comparting of intrinsic reflection and reflection emulation 
is presented in figure [5] 
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Figure 2: Reflection emulation and intrinsic reflection 



In case of emulation we have additional code without essentially new function- 
ality. If we would like integrate a set of projects by such way then we will 
obtain a lot of lumber. Such disadvantage have any chaotic way of integration. 
But chaotic way has top-down design alternative. The top-down paradigm is 
presented in figure[31 

Let us explain meaning of figure [3] First of all we have pure abstract library 
devoted to science and engineering. This library is used by derived libraries 
devoted to different brunches of science and engineering. Libraries can be ab- 
stract. Integration can use other engineering software as third party. Libraries 
are being used by integrated products. Green, pig and purple color are used for 
libraries, third party and integrated products respectively. Italic font is used for 
abstract libraries. 
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Figure 3: Top-down integration of engineering software products 



2 Interoperability issues 

Chaotic way interoperability disadvantage is well known. Chaotic way of inter- 
operability of N products requires development of M^nL) adapters. Adapters 
provide interoperability between objects of any pair of N element set. In figure 
0] (a) we have 8 products. Two sided arrows are adapters. If we use top-down 




a) Chaotic mteiopembuity b) Top-down interoperabuity 



Figure 4: Interoperability 

scheme then we should develop adapters between base object and other objects. 
So we should have N adapters only. The base object is - object in figured] 
(b). 
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3 Category theory as prototype 



Pure abstract library for science and engineering should have very high level of 
abstraction. However this level of of abstraction is already developed in math. 
In mathematics, category theory deals in an abstract way with mathematical 
structures and relationships between them: it abstracts from sets and functions 
to objects linked in diagrams by morphisms or arrows. High level of abstraction 
of Category Theory is explained in [I] . This book contains figure[5]with following 
text. 




Figure 5: Object and arrow 



We did not actually say what a and / are. The point is that they can be anything you like, a 
might be a set with / its identity function. But / might be a number, or a pair of numbers, 
or a banana, or the Eiffel tower, or even Richard Nixon. Likewise for a. 

Let us consider application of this abstraction. Database diagram and control 
systems diagram are presented in figure [6] 
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a) Database diagram (Visual Studio] 
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b) Control system diagram (Simulirtk) 



Figure 6: Database and control system diagrams 



Common features of these diagram are objects and arrows as abstract objects. 
Figure [7] presents database and control systems diagram in single framework 
(There is no one to one map). 

Now we would like to integrate both diagrams. Integrated diagram is presented 
in figure [8] Presented on [8] diagram has following meaning. Database contains 
recorded input signals of control systems. This diagram enable us to define 
response of control system. 
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a) Database diagram 




b) Control systems diagram 



Figure 7: Database and control system diagrams in single framework 
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4 Present day state 



Explanation of idea usefulness is difficult without any implementation. This 
idea shall find further development. However a lot of features are already im- 
plemented. These features are presented in figure O Presented in figure [H 



OPC Foundation components ILonwebservices XNA Game studio ESRI Components 




Figure 9: Hierarchy of assemblies 



scheme reflects only part of implemented features. Source code can be down- 
loaded from http://www.mathframe.com. It is worth to note whole product is 
not a sum of presented in figure O features only. These features can interact and 
interaction make product much more powerful. This thesis will be shown below 
by examples. This hierarchy can inspire some questions. One of reasonable 
question is: "Why we need abstract physical field library?" Now this software 
supports physical fields in 3D space. However this software is declared as univer- 
sal and should support 2D physical fields in future. And 2D fields library shall 
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be inherited from abstract physical field library. Universality had been provided 
from the very beginning of development. Demo applications show prospects of 
further development. 



5 Elementary examples 

Examples of this section are not not related to real science and engineering 
problems. These examples are rather textbook which is very facile for idea 
explanation. 



5.1 Physical fields 

One physical phenomenon acts to another one and there is backward depen- 
dence. This thesis can be exhibited by example presented in figure [TU] We 



//////////////// 



Figure 10: Charged balls 



have two charged balls Ball 1 and Ball 2. Electrostatic force which act to 
Ball 2 depends on relative position Ball 2 with respect to Ball 1. Otherwise 
the relative position as time function depend on motion of Ball 2 and there- 
fore relative position depends on electrostatic force. Common consideration of 
6D motion and physical fields make software much more effective. Now we 
would like simulate motion of Ball 2. Simulation scheme is presented in figure 
[Til Left part of this picture represents electrostatic field of Ball 1. Objects 




Figure 11: Charged balls motion simulation 
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Coord, Field calculator and Field vector provide calculations of next for- 
mula of electrostatic intensity: 



where r is 3D vector of relative position with respect to charge, k is coeffi- 
cient which depends on charge of Ball 1 and system of units. This calcula- 
tion is separated for performance issue. Object Coord calculates -^j, object 
Field calculator calculates components of 3D vector jp^ and Field vector as- 
semblies these components into vector. All these calculation are used by Field 
object. Properties of the Field objet are presented on figure fT2l 
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Figure 12: Properties of Field object 



These properties have following meaning. First of all field parameter is output 
vector of Field vector. Figure [T3l explains meaning of "covariant" term. If 
3D vector is not covariant then its components depend on sensor position only. 
Covariant vector components depend on both orientation and position. Values 
of components are projections of geometric vector to sensor's axes of reference. 
FigurefT3lpresents two orientations of sensor: blue and green. Projections of field 
vector A are different for these different orientations. So we have 3D covariant 
physical field. The Sensor object is a sensor of the field. The Scalars object 
provides components of field vector obtained by the Sensor object. These 
components are used in following motion equations of Ball 2: 
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Figure 13: Covariant field 



where x, y, z are coordinates of Ball 2 and V x , V y , V z are components of its 
velocity. Coefficient a depends of charge of Ball 2 and system of units. These 
motion equations are contained in Motion object. Then coordinates x, y, z 
are used by moved reference frame Sensor's frame. The Sensor object is 
installed on Sensor's frame. So we have backward dependence. Motion and 
therefore current position of Sensor's frame and Sensor depends on current 
field value. Otherwise current field value depends on current position. Let us 
make situation more complicated. New situation is presented on figure 1X41 Here 
we have installed Field on Field's frame and added new object Relative. The 
Relative defines relative motion parameters of Sensor's frame with respect 
to Field's frame. Relative distance is indicated by Chart object. Physical 
picture remains the same. We have added indication only. Now let us install 
Ball 1 on moved platform as it is presented in figure H~5l 
Simulation of this situation is presented on figure [TBI 

Now Field's frame is moved and Motion of ball 1 object defines motion low. 
The time dependence of distance between Ball 1 and Ball 2 has been changed 
by following factors. First of all now Ball 1 which is source of Field is moved. 
Its own motion influence on relative distance. Secondly motion of field source 
influence on field intensity near Ball 2. So absolute motion of Ball 2 will be 
changed. This factor also influence on the distance. These facts are physically 
evident. I would like to exhibit facilities of such approach from point of view 
of software user. The user simply has added motion to Ball 1 and all other 
dependent factors had been automatically taken into account. This facility is 
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Figure 14: Charged balls motion simulation with relative measurements 
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Figure 15: Charged balls with moved ball 1 
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Figure 16: Charged balls motion simulation with moved ball 1 
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impossible without integration. 



5.2 Space aerodynamics and digital image processing 

The key feature of space aerodynamics is that spacecraft interacts with molecules 
which do not collide with each other. Therefore aerodynamic force depends on 
visible area and does not depend on other parameters of spacecraft shape. We 
can use digital image processing for space aerodynamic calculation. Let us con- 
sider it for space aerodynamics. The image of the Mir orbital station is presented 
on figure H71 This photo should be filtered for space aerodynamics usage. This 




Figure 17: The Mir orbital station 



filtering (digital image processing) is presented on figure [T8l The Prototype 
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Figure 18: The Mir orbital station digital image processing 



object contains source image and Result contains filtering result. Other objects 
contain necessary math. This sample presents main advantage of integration. 
Software devoted to space technology only is not effective without digital image 
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processing. Otherwise digital image processing is not effective without advanced 
math. 



5.3 Algebraic topology 

Algebraic topology is not yet everyday tool of engineer or scientist. Now this 
branch of science looks rather as exotic. However exotic tasks are good tests 
for integrability and universality. The "Modern math" has been developed as 
test of prospects of Category Theory approach. In this chapter calculation of 
different topological invariants [2] is considered. Let us consider Klein bottle K 
and real projective space MP 3 (figure [TUl Homology groups of these spaces can 





Klein bottle Projective space 

Figure 19: Klein bottle and projective space 



be calculated as homology of following chain complexes [5] : 

<h s~l ®2 si ^1 



c 2 



Co; 



Chain complex of Klein bottle K and real projective space is presented below: 

z — ^— *• z — ^— ► z — z 

The application representation of it is shown in figure [20l 



Figure 20: Klein bottle chain complex representation 



Chain complex of real projective space M.P 3 is presented below: 

(!) (i.i) 
Z y Z x Z — -A Z 



Chain complexes of both spaces are presented in figure [2T1 
Now we can calculate homology and cohomology groups, homology with coef- 
ficients and other invariants. More details you can find at Category Theory 



project http : / / categorytheory . sourcef orge .net/ 
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Figure 21: Klein bottle chain complex representation 



6 Realistic samples 

Here we consider more complicated samples which are related to real engineer- 
ing problems. This problems include following disciplines: 

- advanced mechanics; 

- processing of signals; 

- statistics; 

- system identification; 

- control theory; 

- celestial navigation; 

- astronomy; 

- geomagnetism; 

- space aerodynamics; 

- digital image processing; 

- virtual reality. All these samples can be downloaded from http : //www. codeproject . com/KB/ architecture/, 
6.1 Advanced mechanics 

Space technology provides good samples of advanced mechanics. Orbital station 
(figure [17]) is a very complicated mechanical object. It is not rigid body. It has 
solar cell panels. These panels are elastic. The station is stabilized by gyros. 
Moreover station configuration is not constant. The "Mechanical aggregate" 
library had been developed for simulation of similar aggregates. The library 
contains aggregate designer. Why aggregate designer? Indeed mechanical equa- 
tions are well known long time ago. But software development for simulation 
of complicated mechanical objects is not quite easy task. Aggregate designer 
make this task much easier. Aggregate designer is integrated into framework. 
This fact enables us provide interpretability of mechanics with physical fields. 
So it is easy to simulate action of magnetic fields on mechanical objects. I will 
consider this task below. Now we would like to create mechanical model of 
spacecraft from models of its modules. Typical spacecraft module is schemat- 
ically presented in figure [22] This module has own coordinates system OXYZ. 
Also it has places of connections. We can connect other modules to this module. 
Behavior of module is defined by following kinematic parameters: 

- radius vector r; 

- velocity V; 

- orientation quaternion Q; 

- angular velocity u>. 
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Figure 22: Spacecraft module 

But module is not rigid in general. And these parameters are not parameters 
of module. These parameters are rather parameters of one point of module. In 
this article we suppose that these parameters are parameters of origin of OXYZ 
coordinates system. Since module is not rigid it has additional degrees of free- 
dom. These degrees of freedom can be interpreted as generalized coordinates 
Qi (i= l,...,n). Instant state of module is defined by following parameters: r, 
Q, qi, q n , V, u, q±, q n . This parameters will be named state variables. 
Parameters V, u>, q\, q n will be called accelerations. Mechanical equations 
define accelerations by state parameters. Accelerations near connection can be 
defined by following way: 

LOi =ei + P Ui V + Q Ui Cj + ^R Uik q h - (1) 

k 

Vi = ai + P Vi V + Q Vi cb + J2 R v t Jk- (2) 

k 

where i is number of connection. Other variables are matrixes which depend 
on state variables. Connection of two modules is presented on figure [23] Both 




Figure 23: Connection of two modules 



modules have equal acceleration near connection. First module acts to second 
one by force F12 and mechanical momentum A/12. Similarly second module acts 
to first one by force F21 and mechanical momentum M21. These forces and 
momentums satisfy following Newton equations: 

F21 — —F21] (3) 
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M 2 i = -M 2X . (4) 
Mechanical equations of module can be represented by the following way: 

V = a + J2 D v i F i + J2 E v i M i ; (5) 

i i 

uj = a + J2 D «>i F i + J2 E »* Mi > ( 6 ) 

i i 

q k = a + Y^ K ki Fi + ^ L kl M,. (7) 

i i 

In these expressions accelerations are independent variables. Fi (Mi) is force 
(mechanical momentum) of i- h connected module. Other vector and matrix 
parameters depend on state variables. Suppose that m - h connection place 
i-h module is connected to n - h connection place j-h module. Then we have 
following evident relations: 





(8) 




(9) 


Fi m — —Fj n , 


(10) 


= Afj n . 


(11) 



Expressions {J) - (fTTj) are linear by accelerations system of equations. This sys- 
tem is fully defined and therefore can be solved. In result we have mechanical 
equations of whole aggregate. It is worth to note that this system is redun- 
dant. If ni (na) is number of freedom degrees of first (second) module then 
system has n\ + ri2 degrees of freedom. However aggregate has n\ + na — 6 
degrees of freedom. The exist software version that avoid this redundancy. But 
I will not describe it in this article. Aggregate designer provides equations of 
full aggregate automatically. Coefficients of ([T]) - (TTTj) equations are inputs of 
aggregate designer. In general these coefficients can depend on states of mod- 
ules. There exist a lot of variants. But we can abstractly define calculation of 
these coefficients and then provide different implementations. Following types 
of mechanical modules are implemented: 

- absolutely rigid body; 

- elastic console; 

- flywheel; 

Equations of rigid body are well known and are not present them here. In con- 
text of aggregate designer rigid body has following properties: 

- mass m; 

- momentum of inertia J; 

- number of connections; 

- positions and orientations of connections. 

These properties can be edited. Elastic console body is a mechanical system of 
infinite degrees of freedom. Usually math model of this object contains finite 
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degrees of freedom with finite set of valuable harmonic oscillations ??. Ev- 
ery harmonic oscillation can be described by following second order system of 
ordinary differential equation: 

Aq + eq + cq = Q. 

Where Q is generalized force, A, e and c are coefficients. Number of harmonics 
and their properties can be edited. 

Flywheels are used in spin stabilization systems of spacecrafts. Following doc- 
uments contain information devoted to spin stabilization systems: 

- http : //www. f reepatentsonline . com/3767139 .html 

- http : //www. aiaa. org/ content . cfm?pageid=406&gTable=mtgpaper&gID=58741 

- http://adsabs.harvard.edu/abs/1966CosRe. . .4. . 173A 

Flywheel is forced by reversible engine (figure [24]) . 




Figure 24: Flywheel 

Otherwise flywheel acts to engine by momentum M x . Since engine is attached 
to spacecraft this momentum is transferred to spacecraft. This momentum is 
used for spacecraft stabilization. Sine flywheel is rotated we have additional gyro 
momentum M gyro . Gyro momentum can be calculated by following expression: 

Mgyro — X UJ. 

where Jf is inertial momentum of flywheel, Q is angular velocity of flywheel and 
u> is angular velocity of engine (and also spacecraft). Total momentum M is 
equal to geometric sum M = M x + M gyro ; Gyro momentum is undesirable fac- 
tor. Stabilization system should require following condition \M gyro \ << \M X \. 
However since engine acts to flywheel value of is being increased by the time. 
Increasing of M gyro compensated by other devices which acts to spacecraft. In 
this article electromagnetic devices will be considered. 

Now we can assembly simulation model of spacecraft presented on figure 1251 
Numbers 1-5 are numbers of connections places of spacecraft. Flywheels at- 
tached to 3, 4, 5 connection places realize angular stabilization of spacecraft 
with respect to axes X, Y, Z. Besides flywheels spacecraft contains three elec- 
tromagnets for stabilization. Mechanical model of spacecraft is presented on 
figure [26] Numbers 1-5 of links are numbers of spacecraft connections. 
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Figure 25: Spacecraft with two elastic consoles and three flywheels 




Figure 26: Simulation model of spacecraft with two elastic consoles and three 
flywheels 
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6.2 Virtual vibration test 



Now we have got simulation model of spacecraft. However this model is not 
facile for development of stabilization system. Control systems theory usually 
uses linearized models. Such model can be obtained by virtual vibration test. 
Note that presented on figure 121)1 is not facile. It contains a lot of big squares. 
Presented software supports compact representation of models. So we can re- 
place model (figure by its compact representation (figured?]) During virtual 




Figure 27: Compact representation of spacecraft simulation model 



vibration test or virtual spacecraft is being forced by momentum which satisfies 
PID control law [I] of momentum has been used: 

M x = M (t) + K-iw + K 2 <p + K 3 J (pdip. 

Where M x is mechanical momentum, is X - projection of spacecraft angular 
velocity, of spacecraft, is rotation angle of spacecraft with respect to X - axis. 
Parameters K\, K 2 , K% are constants. The test purpose is definition of spacecraft 
transformation function. Its definition can be obtained by response on harmonic 
input. The chirp input signal has been used: 

M (t) = C'sm(at + bt 2 ). 

Virtual vibration test scheme is presented on figure [28] The scheme contains 
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Figure 28: Virtual vibration test scheme 



spacecraft's mechanical model and additional math which is necessary for virtual 
vibration test. In result of vibration test we have obtained response that is 
presented on figure |2"!)1 
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Figure 29: Virtual vibration test scheme 
6.3 Nonparametric identification 

Nonparametric identification methods include definition of gain response and 
frequency response. Universal software can easy resolve these tasks. Full expla- 
nation reader can find article devoted processing of signals. Results of digital 
processing is presented in figure [30l 
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Figure 30: Result of digital processing of virtual vibration test signals 



6.4 Parametric identification 

Parametric identification methods include definition of transfer functions. Con- 
trol systems specialists use logarithmic scale for frequency response. This func- 
tion provides clear picture of control object. So we transform functions of previ- 
ous sections to logarithmic scale. Logarithmic transformations of gain response 
a frequency response (figure I30[) are presented on figure 1311 The Y- axis of fre- 
quency response is also logarithmic. Control system specialist can define that 
such charts correspond to following transfer function: 

= k T lS 2 + T 2 s + lT 5 s 2 +T 6S + l 
[S) s T 3 s 2 + T 4 s + 1 T 7 s 2 + T 8 s + 1 ' 

But parameters fc, Ti,...,Ts are unknown. These parameters can be defined by 
nonlinear regression. Regression algorithm is presented in figure [32l Charts in 
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Figure 31: Approximation results 



the left part of figure [3"2"l represents approximated functions (Frequency response, 
sine and cosine 

of phase). Other squares contain necessary math. Approximation result is 
presented in figure |3"T1 

6.5 Celestial navigation 

Any control system requires sensor. We use celestial navigation sensor for space- 
craft control. The sensor enables to define orientation of spacecrafts. There are 
a lot of types of celestial navigation sensors. Here one of possible schemes is 
provided. Suppose that we have equipment that provides celestial images and 
star catalogues. Comparison of image and catalogue enable us to define orien- 
tation of equipment. So we can define orientation of spacecraft. Algorithm of 
this sensor is presented in figure |3"51 First of all let us consider image process- 
ing. Any celestial image contains interfering information. We need filtration 
for its exclusion. Nonlocal digital image processing is being used for this pur- 
pose. This scheme (figure l33|) contains Initial image (Source image) obtained 
by equipment. Little squares in figure l33l provides necessary math. It result we 
have Filtered image (Filtration result). Both images are presented in figure 
IM1 Figure [551 explains image filtration algorithm. If we have 9 closed white pix- 
els then we replace it by one black pixel. Every other pixels are white. Result 
of filtration enables us to obtain X and Y coordinates of black pixels. Then 
this numbers will be compared with star catalogue. The Stat component ex- 
tracts this numbers from Filtered image. Star catalogue is stored in database. 
Necessary information can be extracted by SQL query. Query statement is pre- 
sented below: 
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Figure 33: Celestial navigation algorithm 
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Figure 34: Image nitration in celestial navigation algorithm 




Figure 35: Image filtration in celestial navigation algorithm 
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SELECT RAdeg, DEdeg FROM hip_main WHERE RAdeg > @RAMIN AND RAdeg < ©RAMAX AND DI 
ODEMIN AND DEDeg < @DEMAX AND BTmag > @BTMIN ORDER BY RAdeg 

This statement has following meaning. First of all we consider limited area of 
sky. Declination and right ascension belong to small intervals. Secondly we con- 
sider only susch stars which magnitudes exceed defined constant (in this sample 
the constant is equal to 9). Query result provides following chart presented in 
figure [36] 
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Figure 36: Star catalogue query chart 



We would like compare this chart with filtered image. This operation requires 
a set of math transformations. Essential feature of these transformations is 
Euclidean transformation: 

x' = (x + a) cos ip+ (y + b) sin ip; 

y = — (x + a) sin <p + (y + b) cos ip; 

Parameters a, b, and ip are unknown. Comparison of star catalogue and filtered 
image enable us to define these parameters. Using these parameters we can 
define orientation of spacecraft. 

6.6 Motion of mass center 

Since we consider magnetic method gyro momentum of compensation we should 
have model of Earth's magnetic field. But Earth's magnetic field induction 
depends near spacecraft depends on spacecraft position. So we need consider 
motion of spacecraft as a point. Here I consider usage of inertial reference frame. 
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Figure 37: Inertial reference geometry 



It is more convenient for some tasks. Relation between these reference frames 
is shown in figure [37] 

The OXYZ is inertial frame and OX'Y'Z' is Greenwich one. Greenwich frame 
is rotated with respect to inertial one. We should transform acceleration vector 
g from Greenwich frame to inertial frame. Calculation of g in inertial reference 
frame has two features. Satellite coordinates with respect to OXYZ are not 
equal to coordinates with respect to OX'Y'Z'. Moreover projections of g on 
OXYZ axes of coordinates are different to projections on OX'Y'Z' axes of co- 
ordinates. Declarative approach enables us to resolve both problems at once. 
Usage of covariant fields provides solution of both problems Simulation of linear 
satellite motion is presented in figure |3"81 

Here Earth's Frame component represents Greenwich reference frame. The 
Gravity component is obtained from dynamically linked library. This compo- 
nent calculates Earth's gravity field/ Here this component is used as physical 
field (in Earth's fields). We can say the same about Atmosphere which cal- 
culates parameters of dynamical atmosphere model. Physical field is linked to 
Greenwich reference frame. Gravity field is marked as covariant. These cir- 
cumstances provide solution of both above problems. The Sensor is linked to 
Spacecraft frame. Its orientation coincides with orientation of inertial reference 
frame. Sensor results are used in Motion equations of spacecraft. Otherwise 
Motion equations results are used by Spacecraft frame. So we have math model 
of spacecraft linear motion. 

The picture in figure [38] is not facile for further development since it contains a 
lot of squares. Some of squares should be encapsulated. Encapsulation means 
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that we make invisible some elements in figure [38] Now we would like add model 
of Earth's magnetic field. Result of this operation is presented in figure [39] We 




n 

Figure 39: Encapsulated of spacecraft motion model and magnetic field model 



have full mechanical model. Next step is development of stabilization system. 
Development of stabilization system is based on results which was considered 
in 16.41 Physical principles of control system will be considered in 16.71 

6.7 Final result 

Since this article is more illustration of usefulness of advantages of top-down 
paradigm we leave details of full construction. More details you can find in 
http : //www. codeproject . com/KB/architecture/grandiose2 . aspx. Figure 
l40l shows final result. In brief figure [40] contains two main parts. Left part 
contains full motion model of spacecraft. Full motion model includes: 
-mechanical model of spacecraft; 
-model of Earth's gravity field; 
-model of Earth's atmosphere; 
-model of Earth's magnetic field; 
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-model of Celestial navigation; 

-model of stabilization system of spacecraft. 



Stabilization system could not use flywheels only (See 16. ip . Otherwise there 
is an obstacle for construction of stabilization system which uses electromag- 
nets only. Electromagnetic momentum is always perpendicular to magnetic 
induction B. But stabilization requires all directions of control momcntums. So 
space technology uses double-loop systems which use both flywheels and elec- 
tromagnets. Let us consider both loops of this system. First loop (I name it 
high frequency loop) is presented in figure[5T]We use celestial navigation for def- 



Figure 41: High frequency loop 



inition of orientation and optical gyroscope for definition of angular velocity. In 
result we have following parameters tp x , tp y , <p x , u> x , u> y , lo z . First three param- 
eters are angle deviations with respect to desired axes X, Y, Z. Following three 
parameters are angular velocities with respect to same axes. Control momen- 
tums are obtained by flywheels. We have identified spacecraft angular motion 
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model in 16.41 In accordance to identified model we develop control law. I drop 
details here. Here I note that control law of high frequency loop is designed in 
compliance with control theory. Second loop scheme is presented in figure 321 
This loop purpose is limitation of flywheels' angular velocities. Sensor of this 



loop are tachometers of flywheels. Necessary momentum is provided by electro- 
magnets . It is possible to use different control laws for this loop. But main idea 
of these laws is single In this article I have used following control law. Let be 
Hgyro is total angular momentum. Then momentum of electromagnets is oppo- 
site to H gyro . But electromagnets are not always switched on. If \H gyro \ is too 
small then electromagnets are switched of. Otherwise if angle between Earth's 
magnetic induction H gyro and H gyro is too small then electromagnets cannot 
provide substantionally large momentum that is opposite H gyro . Therefore if 
angle between H gyro and B is too small then electromagnets are also switched 
off. Since action of electromagnets is not continuous I name this loop low fre- 
quency loop. Here I explain how magnetic momentum reduces angular velocities 
of flywheels. Magnetic momentum causes deviation of spacecraft orientation. 
High frequency loop tries to eliminate this deviation by changing of angular 
velocity of flywheels. High frequency loop tries to eliminate this deviation by 
changing of angular velocity of flywheels. Since electromagnetic momentum is 
opposite to H gyro changing of angular velocities is reducing of their values. 
Right part in figure l40l presents spacecraft mission. The mission is astronomical 
observations. Pure mission is shown in figure [43] Spacecraft contains two tele- 
scopes. Telescope 1 is rigidly attached to spacecraft and do not use additional 
frame. Telescope 2 is rotated relatively spacecraft. So Telescope 2 is installed 
on Rotation frame. Otherwise Rotation frame is installed on spacecraft 
body. Special component has been developed for indication of stars. Properties 
of the component are presented in figure 1441 These properties have follow- 
ing meaning. Star coordinates X, Y, Z are equal to FormulaA, FormulaJl 
FormulaJi (of X Y Z component) respectively. Similarly color and size of star 
indication are defined. This sample requires Astronomy Express project and 
star catalogue (Hipparcos and Tycho). 

Both telescopes observes single collection of stars (Stars component). This 
collection is generated by following way. First of all Query component performs 
SQL query of star catalogue. Then other components perform necessary math 




Figure 42: Low frequency loop 
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Figure 43: Spacecraft mission 



Color Size_Foimiia_3 



Color Sze.Foninila_1 



Figure 44: Visualization component properties 
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transformations. Result of these transformations is used by Stars component. 



7 Conclusion 

Typical way of field of activity contains two stages. During first stage there 
is information acquisition and experience without understanding of fundamen- 
tal principles. Many almost equivalent results arc obtained independently by 
different researchers and/or developers. Then fundamental principles become 
well known and chaotic evolution changes to planned development. For exam- 
ple design patterns of GoF [5] had being developed independently After GoF 
these fundamental principles of software development became well known. I 
think that science and engineering software do not reach second stage yet in 
general. There are good fundamental principles of CAD, CAE etc. But they 
are not principles of all science and engineering software. I would like to show 
usefulness of such universal principles. 
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